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We discuss recent work trying to extract the renormalised quark masses and A, the QCD scale, from dynamical 
simulations in lattice gauge theory. 
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1. INTRODUCTION 

Quark masses and Aqcd are both parameters 
which are not accessible to direct observation, un- 
like quantities such as hadron masses which have 
a very clear meaning both in the continuum and 
on the lattice. Extraction of either quantity from 
lattice data inevitably involves both a matching 
of lattice measurements with continuum measure- 
ments and the use of results from perturbation 
theory. 

We can divide the reliance on perturbation the- 
ory into two classes. In the most favourable cases 
all the perturbation theory needed can be per- 
formed in the continuum. In this case the technol- 
ogy of perturbation theory is very well-developed, 
and we can expect to find fairly long perturbative 
series (often to three loops) for the quantities we 
need. Less favourable is the case when lattice 
perturbation theory is needed. The propagators 
and vertices for lattice perturbation theory are 
much more complicated, making high order cal- 
culations very difficult. However we will see that 
some impressive progress is being made with lat- 
tice perturbation theory too. 

The structure of this paper is as follows. First 
we will discuss renormalisation in general, as this 
is crucial to producing quark masses defined in 
a way that can be related to continuum physics. 
Then we will discuss quark masses, dealing sep- 
arately with the lighter quarks (up, down and 
strange) and (in less detail) with the heavier 
quarks (charm and bottom). 

For the light quarks the method we discuss in 
most detail is the use of the pseudoscalar meson 
masses to determine the quark masses. From the 



lowest-order chiral perturbation theory we know 
that a pseudoscalar meson made of quarks of 
flavour a and b has its mass given by 
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which implies that the pseudoscalar meson mass 
(and especially the ir mass) will be very sensitive 
to the quark masses. Eq. also implies that 
the average light quark mass mi = \(m u + m<j) 
will be easier to measure than the separate u and 
d masses, and in this contribution we will only 
discuss this averaged light quark mass. 

In the second part of the paper we will dis- 
cuss attempts to fix the scale parameter of QCD, 
Aqcd, both from measurements carried out at 
"normal" lattice spacings, using gauge configura- 
tions generated for other projects, and from fine 
lattice spacings, using simulations dedicated to 
the purpose of determining A. 

2. RENORMALISATION AND Z FAC- 
TORS 

With very few exceptions, an operator expec- 
tation value measured on the lattice has to be 
renormalised before we can give a number which 
can be compared with experiment. One can think 
of this as a calibration factor. Estimating Z m , 
the renormalisation factor for the quark mass, is 
one of the more difficult steps in the calculation 
of quark masses, so it is appropriate to discuss 
renormalisation in some detail. 

There are a handful of cases where the Z factor 
is not needed (for example, the conserved vector 
current). In other cases one can do a completely 
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non-perturbative calculation. An example of this 
is Zy 1 the renormalisation constant for the local 
vector current, ipj^ip. We know, from the con- 
servation of baryon number and charge, exactly 
what the correct answer should be when the ma- 
trix element of the vector current is measured in a 
hadron. This can be used to calculate Zy, e.g. PQ. 
This calculation does not rely on any perturba- 
tion theory at all, neither on the lattice nor in the 
continuum. 

Usually we are not this fortunate, and we have 
to calibrate our lattice probe by comparing lattice 
results with continuum perturbation theory re- 
sults To use this method one measures quark 
propagators and Greens functions at a range of 
virtualities. Because these are gauge-dependent 
quantities the gauge must be fixed, the usual 
choice is to impose the Landau gauge. The re- 
lation we use to define the renormalisation of an 
operator O is 

A^=^A l S\ (2) 

where Z$ is the wave-function renormalisation 
constant, defined in a similar fashion by compar- 
ing the lattice propagator with the MS propaga- 
tor 1 . 

Where do we get the MS Greens functions 
which we need for this calculation? Here we have 
to rely on continuum perturbation theory. The 
results we need to interpret our lattice results are 
now available up to 3 or 4 loops [3]. Even the 
continuum perturbation series can be slowly con- 
verging in some cases. Tricks to accelerate the se- 
ries' convergence, such as careful choice of scale, 
may be helpful. 

The comparison between lattice and continuum 
results must be performed at a scale \i that sat- 
isfies 

Aqcd < M 2 < l/« 2 - (3) 

1 Traditionally the process is split up into two stages by 
introducing an intermediate renormalisation scheme such 
as RI or RI'. A renormalisation factor taking us from 
lattice to RI' is followed by a conversion factor to change 
RI' to MS. I think the nature of the calculation is easier 
to follow if we combine both stages, and convert directly 
from lattice to MS. 



The lower limit arises because when is too 
small we do not know the true value of Aq S , 
since our perturbative series do not converge fast 
enough, and because the Greens function may 
have large non-perturbative contributions, for ex- 
ample contributions associated with the sponta- 
neous breaking of chiral symmetry and the ex- 
istence of a quark condensate. The upper limit 
arises because the lattice Greens function will 
have discretisation errors which become impor- 
tant when a 2 fi 2 ~ 1. The method of [5| relies on 
the existence of a plateau region between these 
two limits, and in practice this may turn out to 
be disappointingly narrow. 

There is little that can be done about the lower 
limit, what happens here is a real physical effect, 
present both on the lattice and in the continuum. 
On the other hand lattice perturbation theory can 
be used to greatly reduce the a 2 /! 2 errors, extend- 
ing the useful plateau region outwards. 

Normally lattice perturbation theory results 
are quoted for small external momenta a 2 /i 2 <C 1 
but with a little extra effort the Greens functions 
can be found for any external momentum. Sub- 
tracting off the one-loop contribution to the dis- 
cretisation errors allows us to reduce errors of this 
type from 0{g 2 a 2 ^ 2 ) to 0(g 4 a 2 /j, 2 ). The success 
of this procedure can be seen in Fig. ^ which 
shows a much better plateau after subtraction of 
one-loop lattice artefacts. 



2.1. Singlet Contribution 

We have discussed how to measure the renor- 
malisation constants for quark bilinear operators, 
but how can we measure the mass renormalisation 
factor Z m l The usual method is to look for an 
identity relating Z m to Z$, the renormalisation 
factor for the scalar operator. 

The usual way to derive this is by looking at the 
divergence of a flavour non-singlet vector current 
J°. In the continuum this current is ^yj "^, 
where r is a traceless flavour matrix. The lattice 
form depends on the fermion action, for Wilson 
or clover fermions it is 
i^(x + /t)[7 M +r]r a ^(a;) 

+ ^(»bv ~ r]U^(x)T a ip(x + A). 
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Figure 1. The renormalisation factor for the 
axial current, for quenched clover fermions at 
(3 = 6.20. The open squares show the raw data, 
the black squares show the result after subtraction 
of the one-loop discretisation errors. The subtrac- 
tion procedure increases the plateau region. 

The Ward identity 2J reads 

(dfj.jfy = (V> [M,r a ] V) + contact terms, (4) 

where M is the quark mass matrix. This iden- 
tity links the flavour non-singlet scalar operator 
and the flavour non-singlet part of the mass ma- 
trix. Eq. (0J is an exact lattice identity if we use 
the conserved vector current, but because it only 
involves the non-singlet part of the mass matrix 
(i.e. mass differences), it still permits an additive 
renormalisation of the singlet part of the mass 
matrix (which does of course occur in Wilson or 
clover fermions). 

Requiring that eq. Q still hold after renormal- 
isation tells us that Z% S Z^ S = 1. This tells us 
how to renormalise non-singlet quantities, such 
as mass differences, but it doesn't tell us how the 
singlet part of M (which commutes with r) renor- 
malises. 

In MS we are used to having a difference be- 
tween singlet and non-singlet renormalisation fac- 
tors for operators with an odd number of gamma 
matrices. Because clover fermions have no ex- 
act chiral symmetry, there can be a difference be- 
tween singlet and non-singlet Z for operators with 
an even number of gamma matrices too. This can 
be seen in Fig. [21 In the continuum the "bubble" 



diagram would be the product of 5 gamma matri- 
ces (3 from propagators, 2 from qqg vertices) so 
its trace would be zero in the massless case. For 
Wilson or clover fermions this no longer holds, 
and the bubble will give a non-zero contribution. 





Figure 2. An example (above) of a Feynman dia- 
gram present in both singlet and non-singlet cases, 
and (below) a diagram which only contributes to 
the singlet Z. The black point shows where the 
ipip operator is inserted. 

The singlet vector current does not give us 
any useful information about the mass matrix, 
but we can find some relationships by considering 
partially quenched QCD, with valence quarks al- 
lowed to have a different mass from the sea quark 
mass. If we change the valence quark mass while 
keeping the sea quark mass fixed, the derivative 
of the quark propagator S gives the non-singlet 
scalar three-point function 
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-S(p) 
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(5) 



However if we change valence and sea quark 
masses together we get an identity for the singlet 
scalar three-point function 



d 



d 



2k sc , 



S(p) = G^(p). (6) 



From these identities we conclude that 
Z^ Z^ — i (which we already knew from the 
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Vector Ward Identity) and Z^Z— = 1 which 
is new. We can find Z^ s and thus Z^ s by the 
non-perturbative methods discussed earlier. The 
identities here also give us a way of finding Z^. 
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Figure 3. Dynamical and partially quenched pseu- 
doscalar meson masses for clover fermions at 
P = 5.29. Data from JJ%. 

In Fig. |21 we show partially quenched (open 
points) and full dynamical (filled points) data for 
the pseudoscalar meson masses. We can immedi- 
ately see that there is a different critical k if we 
follow the solid line connecting the full dynami- 
cal points (the points with valence and sea quark 
masses equal) than if we follow one of the dashed 
lines (sea quark mass held fixed, valence quark 
mass varied). This is because clover fermions 
have no exact chiral symmetry, and so an additive 
quark mass renormalisation is possible. The dif- 
ference in critical us shows that this additive term 
depends quite strongly on the sea quark mass. 

In the continuum the situation would look quite 
different — because of chiral symmetry both the 
full and partially quenched curves would both 
have to pass through the origin, at zero valence 
quark mass. This would also apply for a lattice 
fcrmion formulation (such as overlap fermions) 
which have an exact chiral symmetry. After 
renormalisation the lattice results ought to show 
the same structure, with both full dynamical and 



partially quenched mpg vanishing at the same 
place. The only way to arrange this is to use 
different renormalisation factors for the partially 
quenched and full QCD quark masses, see Fig.^J 
just as we expected from considering the identi- 
ties eqs. (0 and ©. 
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Figure 4. Dynamical and partially quenched 
pseudoscalar meson masses before and after 
renormalisation of the quark mass. 



We define the bare sea and valence quark 
masses by 

m=- (7) 

2/c 2/^ c 
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with k c defined as the critical k for QCD with 
equal valence and sea quark masses (see Fig. 0}. 
The renormalised masses are then defined by 



sea 



Z^msea , (8) 
( m val ~ m sea) + ^m sea . 

The partially quenched pseudoscalar mesons are 
massless at the point where m^ al — 0, which we 
call n v c al . We can use this to find the ratio of the 
two Z's, 

Z~~ lYlsen. TH"m 



Z NS 



(9) 
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Z ratios calculated this way in [S] are shown in 
Fig. [5] The ratio only depends weakly on the 
quark mass, but it depends very strongly on 0. 
The ratio drops rapidly as lattice spacing de- 
creases. 



flavours of unimproved Wilson fermions. In this 
paper they initially found a value for the ratio 
m s /mi « 52, which is a long way from the value 
found in quenched calculations, and also far from 
the predictions of chiral perturbation theory 



2m 2 K 



mi 



25 . 



(10) 



Normally one would expect the mass renormalisa- 
tion factor to cancel out in this ratio, but because 
the singlet and non-singlet Z differ this is not en- 
tirely true. 

In Fig. E| we show a pion line, the line giv- 
ing the pseudoscalar mass when valence and sea 
quark masses are equal, and the kaon line, show- 
ing the result when sea quarks and one valence 
quark (representing the u or d in the K) are kept 
at a fixed mass, while the other valence quark 
(representing the strange quark) is varied. The 
slopes of these two lines are different, just as seen 
in the partially quenched case. 
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Figure 5. The ratio Z^/Z^ for dynamical clover 
fermions (3 values run from 5.20 (highest 

line) to 5.40 (lowest line). 



Figure 6. A schematic diagram illustrating the 
subtleties of measuring the strange/light quark 
mass ratio from Wilson fermions. 



One way to illustrate the importance of the ra- 
tio Z^/Zm is to go back to one of the pioneering 
dynamical calculations 6 , carried out with two 



We determine the k corresponding to the light 
quark mass by seeing where the it line crosses the 
mass-squared of the real physical pion, and the 
strange quark k by seeing where the kaon line 
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crosses the mass-squared of the physical K. If 
we now use these two k values to calculate the 
ratio of the bare quark masses, m b s are jm\ are we 
get a ratio much larger than from eq. (|1U|> . be- 
cause the K line in Fig. H3 is less steep than the 
7r line. However, if we renormalise the masses ac- 
cording to eq. (BJ, then the slopes of the kaon 
and pion lines will become approximately equal, 
as in the lower panel of Fig. 21 and now the ratio 
of the renormalised quark masses will be close to 
eq. JHJ. 



a, so one will not be able to completely remove 
this effect by power-law extrapolations in a. 

The argument we give here, using different Z 
factors for the singlet and non-singlet parts of the 
mass matrix is rather similar to the discussion 
given at Lattice '97 in 

2.2. Quark mass results 

The preliminaries are now over, and it is time 
to discuss some of the recent quark mass results 
from dynamical simulations. 
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Figure 7. A figure from J7y showing the 
ence in light quark masses found from partially 
quenched and fully dynamical ("direkt") data. 



In later works the SESAM collaboration 
used a different strategy, measuring data from 
both partially quenched and fully dynamical 
mesons, and then making a joint extrapolation 
to the continuum limit, Fig. The same Z was 
used in both cases. This gives a result more in 
line with what we do here, but it is not exactly 
the same. Although we have seen in Fig. that 
the ratio Z^/Z^ s decreases as a decreases, there 
will presumably be a perturbative contribution at 
the two-loop level, which means that the Z ratio 
will asymptotically approach the value 1 rather 
slowly, like 1 + 0(g i ), rather than like a power of 
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Figure 8. A comparison of quenched (above) and 
dynamical (below) estimates for the strange quark 
mass 1$). 



In Fig. |3 we show an estimate of the strange 
quark mass calculated using the flavour-singlet 
renormalisation constant calculated from the 
comparison of dynamical and partially quenched 
mesons. The result is rather similar to that found 
in the quenched case. If the non-singlet renormal- 
isation constant had been used the a-dependence 
would be much stronger, and the data points 
would have been lower. 

An alternative mass definition is based on the 
continuum Ward identity for an axial current 
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with flavour index a 0], 

(df.A^) = (^{m,T°} 75^) + contact terms. (11) 

In clover fcrmions the axial current has to be im- 
proved by adding irrelevant operators, because of 
the lack of a true chiral symmetry. As can be seen 
from the Ward identity, this axial Ward identity 
(AWI) mass has to be renormalised by the factor 



N^2+1 



N^2 



Za AWI 
——TO 



(12) 



The final renormalised result ought to be the 
same in either case, but we will see that this is 
not yet the case, there are still fairly large dif- 
ferences between different groups working on this 
problem. 

The AWI quark mass definition has been used 
by the CP-PACS and JLQCD Collaborations to 
study the quark masses in clover QCD, both for 
two- flavour [S] and (2+1) flavour simulations |10j . 
The Z ratio needed to renormalise the bare lat- 
tice masses are calculated from tadpole-improved 
one-loop lattice perturbation theory. Their re- 
sults are summarised in Fig. [5] In the quenched 
case the result depends quite heavily on whether 
the K or (f> meson is used to determine the strange 
quark mass. This is a sign that the quenched the- 
ory differs from the real world. Fortunately this 
discrepancy is greatly reduced when sea quarks 
are switched on, as can be seen in the panels of 
Fig. showing the s quark mass from dynamical 
QCD. 

The SPQcdR Collaboration |TJ have measured 
quark masses in an Nf = 2 simulation using the 
Wilson quark action at a single (3 value, (3 — 5.8. 
They use both the VWI and AWI definitions 
of quark mass, and find compatible results from 
both methods. The rcnormalisation constants are 
found non-perturbatively using the method of [2j ■ 
Their preliminary values are included in Table ^ 
and are broadly compatible with the results of . 

In |12| the quark masses are measured in a 
dynamical simulation with (2+1) flavours of sea 
quark. The use of staggered quarks speeds up the 
calculation considerably, allowing the light u/d 
sea quarks to be simulated all the way down to 
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Figure 9. The u and d quark masses (above) and 
s quark mass (below) at a ~ 0.1 fm, from 



masses ~ to s /8, a considerable advance on what 
is possible with other quark formulations. The 
staggered quarks have been improved to suppress 
'taste' violation and other discretisation errors. 
This also has the effect of making perturbation 
theory better behaved than it is for unimproved 
staggered fermions. The Z factors used in [T2] 
are calculated in one-loop tadpole improved per- 
turbation theory. The estimated error due to the 
unknown higher-order terms in perturbation the- 
ory is 9%. This study finds that the effects of 
including sea quarks are dramatic, and that the 
final physical masses of the quarks are a lot lighter 
in the dynamical simulation than they were in a 
quenched calculation. 

Fig. EH shows determinations of the strange 
quark mass compiled by the Particle Data 
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mi /MeV 


m s /MeV 


sea quarks 


Reference 


3.05(6) 


80.4(1.9) 


(2+1) 


CP-PACS and JLQCD [TD] 


2.8(0)(1)(3)(0) 


76(0)(3)(7)(0) 


(2+1) 


HPQCD, MILC and UKQCD [E] 


4.7(2)(3) 


119(5)(8) 


Nf = 2 


QCDSF/UKQCD 


4.8(5) 


111(6) 


N f = 2 


SPQcdR Collaboration [TT] 


4.5(17) 




Nf = 2 


SESAM/TxL 



Table 1 

Masses of the light, mi — \(m u + m^), and strange quarks in dynamical lattice QCD, as determined by 
various collaborations. The masses are given in the MS scheme at a scale of 2 GeV. 
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Figure 10. Determinations of the strange quark 
mass from the Particle Data Group 



to have expressions resumming all (am c ) n terms. 
Suggestions on how to do this (known as FNAL 
masses) were made in |14|. As the lattice spac- 
ing becomes smaller all definitions should agree. 
This is examined in |15| . where the charm mass 
determined from the axial and vector Ward iden- 
tity definitions, and the mass determined from 
the FNAL method mi and 7712 are compared in 
the quenched case (where the a range can be 
made large). At present lattice spacings the mass 
from the Axial Ward Identity (m^) looks as if it 
might extrapolate to a lower value than the oth- 
ers, though there is still room for it to move up 
again, see Fig. ^] 

This reference also looks at the effect of find- 
ing m c on dynamical clover configurations, the 
open symbols on the right-hand side of Fig. ITTI 
There seems to be very little difference between 
the quenched and dynamical results, though of 
course the present results are at a single rather 
coarse lattice spacing, and the sea quarks in the 
simulation are rather heavy. 



Group Comparing this with Table ^ we 

see that most of the dynamical lattice determi- 
nations look low compared with estimates from 
traditional methods. 

2.3. The charm quark mass 

The charm quark mass has its own special dif- 
ficulties because m c is comparable with the in- 
verse lattice spacing a -1 for the lattice spacings 
we have to work with at present. This means 
that lattice artefacts proportional to high powers 
of am c can be important, so it would be useful 



2.4. b quark mass 

A nice example of the interplay between pertur- 
bation theory and physical result is the b quark 
mass calculation of Gimenez, Giusti, Rapuano 
and Martinelli 16 converting lattice data into a 
value for mj, in MS, using stochastic perturbation 
theory results from |17| . 

This requires perturbation theory in two places, 
the conversion from the lattice mass parameter of 
heavy quark effective theory (HQET) to the pole 
mass requires the subtraction of a residual mass 
dm which is proportional to 1 /a. This subtraction 
can be calculated in lattice perturbation theory. 
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Figure 11. Different determinations of the charm 
quark mass \15\/ . The closed symbols show the 
effect of using different mass definitions in the 
quenched case. All definitions should lead to the 
same continuum value. The open symbols at the 
right show the charmed quark mass determined on 
Nf = 2 dynamical configurations from UKQCD. 



After that, the conversion from pole mass to MS 
mass is done using results from continuum per- 
turbation theory. Since the b quark pole mass is 
just used as an intermediate step the calculation 
taken as a whole is presumably legitimate even 
though confinement means that there is no pole 
in the 6-quark propagator. 

The two-loop subtraction for the residual mass 
was calculated using traditional lattice perturba- 
tion theory methods in ^Hl • The three-loop result 
was calculated using the novel method of stochas- 
tic perturbation theory, which has the major 
advantage that the difficulty of calculating new 
terms in the series does not increase as rapidly 
as it does in conventional calculations. Stochas- 
tic perturbation theory thus promises to over- 
take conventional techniques in the future. It 
does have some disadvantages, being a stochastic 



method, the coefficients calculated do have a sta- 
tistical error, and can not be given to the number 
of decimal places that we are used to from con- 
ventional perturbation theory, but they are usu- 
ally quite good enough for comparison with real 
lattice data. 

The dependence of the final b quark mass on 
the order of lattice perturbation theory used is 
illustrated in Table |21 Knowing the three-loop 
coefficient helps to significantly reduce the sys- 
tematic uncertainty in the final estimate, which 
is 

m b (m b )\ unq = (4.21 ± 0.03 ± 0.05 ± 0.04) GeV. 

3. DETERMINING A 

An important (but difficult) project for the lat- 
tice is to link the low-energy, non-perturbative 
quantities we normally measure (such as hadron 
masses) with the high-energy world accessible in 
accelerators, where perturbative QCD works well. 
The key to bridging this gap is to find a value for 
A, the scale parameter of QCD. 

The strategy for determining A has two nec- 
essary parts. Firstly some dimensionful quantity 
has to be measured to set the scale. A common 
choice is to look at the static potential. This can 
be measured fairly well. A good way to fix a 
length scale is to find the distance at which 

r 2 F(r) = -r 2 ^-V(r) (13) 
ar 

has a particular value 19 , the usual value chosen 
is 1.65, which gives a distance scale ro ~ 0.5/m. 
Another way to find a length scale from the po- 
tential, now less popular, would be to use the 
string tension. 

Secondly, we need to extract the MS coupling 
at a known scale, which can then be converted 
into a value for hjjg- One possibility is to use 
lattice perturbation theory to relate the lattice 
coupling g 2 = 6//3 to g^. 

Relating the A parameters of two schemes is 
a finite calculation. This is an asymptotic state- 
ment, and can be answered completely by a one- 
loop calculation 20 . 

However, relating the coupling constants at a 
finite scale is no longer so simple. Outside the 
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m b {m b ) 


Dependence on higher orders (unquenched data) 


Order 


5m 


pole 


m b {m b ) / m p b ole 


m b (m b ) 


Am b 


LO 





3770 


1.000 


3770 




NLO 


978.79 


4748.79 


0.92276 


4382 


612 


NNLO 


10.53 


4759.32 


0.89060 


4239 


143 


NNNLO 


74.55 


4833.87 


0.86883 


4200 


39 



Table 2 

A table illustrating the way in which the estimate for the bottom quark mass m b {m b ) depends on the 
order of the lattice perturbation theory used. The quark mass results are preliminary values from ]16f . 
using perturbative calculations from \17Jj and ]lttf . 



asymptotic regime we need calculations of the /3- 
function of both schemes to as high an order as 
possible, this can then be translated into a rela- 
tion between the couplings. Currently we know 
the lattice (3 function to three loops |21| . 

As often happens with lattice perturbation the- 
ory we find that the series converge rather poorly. 
Fortunately we have several methods to accel- 
erate the convergence. One explanation for the 
large coefficients in lattice perturbation theory 
is that they come from gluon "tadpoles" due 
to the higher order vertices present on the lat- 
tice [22] • A way to take these into account is to 
re-express the series in terms of a "boosted" cou- 
pling, «?2 = g 2 /UP- 
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Using the quenched case as an example we can see 
in the equation above that the coefficients in the 
series for the MS coupling are reduced when we 
change the expansion parameter from g 2 to g 2 ^. 
We can improve the series further by choosing a 
natural scale If, instead of g|^-(l/a) we cal- 
culate (/^j.(2.63/a), the term vanishes, and 
the first correction is 0(g 2 ^) with a rather small 
coefficient. This gives us grounds to hope that 
the unknown higher-order coefficients will not se- 
riously change the estimate of g^r- 

Once we have values for g^r we can use the 



known (3 function to convert these into a value 
for A M . The results in the quenched case |24j . 
based on data from [35] and are shown in 
Fig.[H 
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Figure 12. A MS calculated for quenched QCD 
for (3 values in the range 5.95 to 6.92 \24\ - The 
ro values for high (i are from I2fi\l . The data ex- 
trapolates very smoothly towards a value A^^ =0 = 
242(11) (10)MeV. 

Some extrapolation in a is needed, this seems 
not to present any great problem, the data look 
surprisingly linear in a , and the smallest a values 
are close to the c onti nuum limit. The result in 
physical units is A^ = 242(11) (10)MeV. The 
first error is a statistical error, the second error is 
an attempt to estimate the possible effects of the 
unknown higher order terms in the (3 functions. 

One can try a similar calculation for the dy- 
namical case, using the clover action configura- 
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tions produced by the UKQCD/QCDSF Collab- 
orations. The lattice [3 function for clover quarks 
is available In the dynamical case we need 
a double extrapolation, both in lattice spacing 
and in sea quark mass. Both these extrapola- 
tions are over a rather large distance. A simple 
phenomenological extrapolation formula is 

r A MS = const + B(a/r ) 2 + Cam q + Dr m q . 

(15) 

The formula fits the data fairly well, as can be 
seen in Fig. ^] 
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Figure 13. Extrapolation of A MS for dy- 
namical quarks. The data are measured on 
UKQCD/QCDSF configurations with dynamical 
clover fermions \24\ - The upper part of the figure 
shows the mass dependence after the data has been 
extrapolated to the continuum limit. The lower 
half shows the a dependence of data after extrap- 
olating to the chiral limit. 



The extrapolated value hasn't changed very 
much due to the jump to finer lattice spacings 
in the past year [23], the new value is A = 
222(4) (28)MeV at N f = 2. 



Figure 14. A comparison of continuum extrapola- 
tions ofroA MS \27Jj made for various fermion ac- 
tions. Wilson fermions are extrapolated linearly 
in a, while clover and staggered are extrapolated 
linearly in a 2 . 

One way to estimate the systematic error in 
the final result is to repeat the calculation for 
different formulations of dynamical fermions, and 
see how much the results scatter. A very complete 
study was carried out in [57]. As can be seen 
in Fig. [21 the results do still depend quite a lot 
on the action and extrapolation procedure. The 
dependence on the action could be due in part 
to higher order terms in the perturbation series, 
which will of course depend on the fermions used. 

A very impressive attack on finding A was re- 
ported in this year's Lattice Conference [2Hj, ex- 
tending work done in 29 with Asqtad improved 
staggered fermions. Extensive three-loop calcula- 
tions of Wilson loops and the static potential have 
been carried out. Some idea of the work involved 
can be gained by looking at the number of Feyn- 
man diagrams that have to be considered, Fig. 1 151 
Wilson loops of different sizes are dominated by 
gluons of different virtualities, with small loops 
corresponding to large Q 2 . Thus by calculating 
a s from loops of different sizes one can see how 
a s runs even from data taken at a single value of 
a. This is illustrated in Fig. 1161 where one can 
see how well the running of a agrees with the [3 
function predictions. 
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Figure 15. The gluon and ghost Feynman dia- 
grams used in ]2Sf to calculate three-loop results 
for Wilson loops and the static potential. A simi- 
lar number of fermion loop diagrams is needed to 
calculate the unguenched case. 



In Fig. El we see the lattice results extrapolated 
out to the scale mz, where they are compared 
with the traditional Particle Data Group value. 
The agreement looks impressively good. 

Is there any way of diminishing the dependence 
on high-order perturbation theory? The obvious 
way to reduce our sensitivity to high order terms 
in the perturbation theory is to work at weaker 
values of the coupling — in QCD this means 
working at high energy or short distance. This re- 



Running of alpha 




q* in GeV 



Figure 16. The running of a s as determined 
by \2ty . a s is defined in the potential scheme. 



quires a very fine lattice spacing, a, which in turn 
forces us to consider simulations in small physi- 
cal volumes. The small volumes involved in this 
method mean that it isn't possible to reuse config- 
urations generated for other projects on hadronic 
physics, the calculations require many dedicated 
runs at very large /3 values. 

The ALPHA Collaboration have been pursu- 
ing this method for many years |3(J| . using the 
same methods that they applied in the quenched 
case A recent plot of their results is shown 
in Fig.^J By repeatedly halving the lattice spac- 
ing (and the physical size L of their system), 
they make measurements of a s over a large range 
of scales (more than two orders-of- magnitude). 
These are compared with the running expected 
from the three-loop (3 function. The data reach 
a s values ~ 0.1, where we are very confident that 
perturbation theory will work well. As can be 
seen, the agreement with the perturbative f3 func- 
tion remains close, except perhaps at the very 
lowest scale. The fact that perturbation theory 
works so well over such a large range offers some 
comfort to the groups that base their determina- 
tion of A on lattices with "normal" lattice spac- 
ings, who have to rely on data with /i/A around 
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Lattice Results Compared With PDG 
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Figure 17. a s (mz) estimated from Wilson loops 
of various sizes, and from the short distance po- 
tential p%y . 



Figure 18. The running of a s in two-flavour dy- 
namical QCD as determined by the ALPHA Col- 
laboration '301 



10 or 20, but clearly the precision that will be 
achieved from data at fi/A ~ 1000 should be 
much greater. 

This work has not yet given a final value for 
A in physical units, preliminary results give the 
value ]n(AL max ) = -1.34(7) where a s (L max ) = 
0.372, but the final step is still needed, a con- 
version from the length scale they use, L max , to 
physical units. 

4. CONCLUSIONS 

The determination of quark masses is the more 
mature of the two fields we have looked at. We 
see that agreement has not yet been reached, and 
there are still differences between the measure- 
ments of the different groups. It is not yet clear 
(to me at least) where the origin of these differ- 
ences lies, as there are many differences between 
the calculations, particularly in the ways the Z 
factors are calculated and in the choice of mass 
definition. 



One area in which we are making progress, and 
can expect more advances in the coming years, is 
bringing the sea quark masses down towards more 
realistic values. When simulations are carried out 
with large sea quark masses it can be difficult to 
find much difference from quenched calculations, 
and the extrapolation to the physical region can 
be difficult. Now, with more powerful machines 
and new fermion actions, we are seeing simula- 
tions carried out much closer to the physical re- 
gion, and we can expect more soon. 

A topic which we have not covered here is 
the calculation of quark masses using overlap 
fermions, since dynamical overlap is not yet as ad- 
vanced as dynamical calculations with the older 
fermion formulations. 

In both the calculation of masses and of A 
we have a need for higher orders of perturba- 
tion theory. Advances in lattice simulations are 
making lattice perturbation theory more difficult, 
because more complicated actions lead to more 
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complicated Feynman rules, and because there is 
now far greater diversity in the choice of fermion 
and gauge actions used. At Lattice '04 we have 
seen that diagrammatic techniques are rising to 
the challenge, we have also seen that stochas- 
tic perturbation theory is delivering results which 
are useful to those extracting physics from lattice 
simulations. 
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